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Abstract: We consider a mathematical model for the interactions of an elastic body fully 
immersed in a viscous, incompressible fluid. The corresponding composite PDE system com- 
prises a linearized Navier-Stokes system and a dynamic system of elasticity; the coupling takes 
place on the interface between the two regions occupied by the fluid and the solid, respectively. 
We specifically study the regularity of boundary traces (on the interface) for the fluid velocity 
field. The obtained trace regularity theory for the fluid component of the system — of inter- 
est in its own right — establishes, in addition, solvability of the associated optimal (quadratic) 
control problems on a finite time interval, along with well-posedness of the corresponding op- 
erator Differential Riccati equations. These results complement the recent advances in the 
PDE analysis and control of the Stokes-Lame system. 
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1 Introduction 

The description of the interactions between a fluid flow and an elastic structure via 
an interface is a classical topic in continuum mechanics. Motivated by an impressive 
variety of fields of application, such as, e.g., cell biology, biomedical engineering, ul- 
trasound imaging, in the past decade fluid-structure interaction models have received 
a renewed attention. Numerical and experimental studies have pursued the derivation 
of appropriate PDE models for specific physical problems, along with the understand- 
ing of the corresponding dynamics; see, e.g., [26] and [25]. On the other hand, basic 
mathematical questions such as well-posedness (in a natural functional setting) of the 
systems of coupled partial differential equations (PDE) which result from the modeling 
of these interactions, have been explored and answered only recently. 

Let us recall that a major technical issue which arises in the PDE analysis of 
fluid-solid interaction models is the apparent discrepancy between the regularity of 



1 The research of the first author is supported by the Italian MIUR, within the project 2007WE- 
CYEA ("Metodi di viscosita, metrici e di teoria del controllo in equazioni alle derivate parziali non 
lineari"). The research of the second author is supported by the NSF grant DMS-060666882 and by 
the AFOSR grant FA9550-09- 1-045. 



1 



boundary traces for the fluid and the elastic (hyperbolic) components on the interface, 
with the latter not defined a priori by means of standard trace theory. Different solu- 
tions have been proposed to cope with the mismatch between parabolic and hyperbolic 
regularity. For instance, in [9] a regularizing term is included in the elastic equation, 
thereby ensuring the soughtafter boundary regularity, as a consequence of the im- 
proved interior regularity; the works [13] and [15] consider instead very smooth data, 
thus obtaining only local (in time) existence. In contrast with previous work, the study 
of [6j achieves — at least in the case of static interface — well-posedness of an established 
nonlinear fluid-structure interaction model in the natural energy space, along with the 
exceptional yet intrinsic regularity of the boundary traces of the (elastic) hyperbolic 
component of the system on the interface. (To the authors' best knowledge, a similar 
well-posedness result for the more challenging case of moving interface is, at present, 
lacking.) 

The composite PDE system under consideration, which is problem f)2. 1 1) in the next 
section, is a linearized version of the nonlinear system studied in [HE]; from a modeling 
point of view, the linear evolution is consistent with the case of very slow motion of the 
fluid. The boundary value problem comprises a linearized Navier-Stokes system (often 
referred to as Oseen equations) for the velocity field u of the fluid flow and a dynamic 
system of elasticity for the displacement w of the solid. The interaction between the 
two means takes place at an interface between the two regions occupied by the fluid 
and the solid, respectively. We consider the case in which the motion of the solid is 
entirely due to infinitesimal displacements, and therefore the fluid-solid interface, later 
denoted by T s , is stationary. (Yet, although the displacement u is small, the velocity 
ut is not, resulting in the interface condition u = Wf.) While the PDE problem (|2.1j) 
is present in the literature since the late sixties (see [23J ) , its well-posedness was not 
shown until [6], whose theory encompasses the present model. In the case of the 
(linear) Stokes-Lame system, the contemporary work of [3j HI [5] provides not only 
well-posedness, but also a careful stability analysis. A very nice overview of recent 
(theoretical as well as numerical) contributions to this topic is found in [TJ], whose 
introduction also provides insight into modeling aspects; see the references therein. 

The present study focuses on the regularity of boundary traces for the velocity 
field of the fluid flow u (and ut) on the interface T s . Although the obtained regularity 
properties are of interest in their own right, they are also critical in the study of the 
optimal boundary control problems (with quadratic functionals) associated with the 
fluid-solid interaction (|2.1j) , in the presence of active controls on the interface T s , thus 
constituting its main motivation as well as application. 

It was shown in [21] that the fluid-structure interaction model under consideration 
yields a suitable ('singular') estimate for the corresponding abstract evolution, which 
plays a major role in the study of the associated optimal control problems. Indeed, in 
light of the regularity result obtained in [21^ Theorem 5.1], the optimal control theory 
in [17\ \19\ [20] applies, ensuring a feedback control law with a bounded gain operator 
(defined on the state space), along with well-posed Riccati equations. This, however, 
holds true assuming that the functional penalizes the quadratic energy of the solid 
component in a weaker topology than the actual energy level (this is a consequence of 
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the fact that the singular estimate holds in the norm of a function space TL~ a which 
is larger than the energy space). 

A first attempt to remove the constraint on the structure of the quadratic cost 
functional — which corresponds to a constraint on the observation operator — has been 
successfully carried out in [T2] in the case of the simpler Stokes-Lame system. The PDE 
analysis performed in [12] shows indeed that the singular estimate in [211 Theorem 5.1] 
does not hold for a = and hence that the optimal boundary control theory of 
[171 [20] is inapplicable to the case of natural cost functionals such as the integral of 
the quadratic energy of the system. Nonetheless, [12] provides novel regularity results 
which ensure solvability of the associated optimal control problems, with well-posed 
Riccati equations, according to the distinct (optimal control) theory contained in [2]. 

In the present paper we show that the complex of boundary regularity results estab- 
lished in [12J for the uncontrolled Stokes-Lame system may be extended to more general 
fluid-structure interaction models comprising a linearized Navier-Stokes (rather than 
the Stokes) system for a slow, viscous, incompressible fluid (see (I2.ip below). As we 
shall see, the regularity of boundary traces for the fluid component of the system will 
follow combining the exceptional boundary regularity of the hyperbolic component 
of the system (c/. Lemma 12. TH . the regularizing effect of the parabolic component, 
and interpolation techniques. The same factors are behind the proof of the specific 
regularity estimates recently derived for diverse composite PDE systems, including 
thermoelastic systems ([HJQ]) an d acoustic-structure interaction models (|10j). 

An outline of the paper follows below. In Section [2] we introduce the coupled PDE 
system, along with the natural functional setting where well-posedness holds true, 
in both the variational and semigroup sense. The statement of the present work's 
main contribution, namely Theorem I2.6[ which collects several boundary regularity 
results for the fluid velocity field of problem (|2.ip . is given here. Section [2.41 contains 
some preliminary trace regularity results which are crucially instrumental to the proof 
of Theorem 12.61 Section [3] is entirely devoted to the proof of all the assertions of 
Theorem 12.61 Application of the obtained regularity theory to the associated optimal 
control problems is discussed in Section [H which concludes the paper. 

2 The mathematical model, statement of the main result 

In this section we introduce the PDE system under investigation. We preliminarly 
recall from [6] the well-posedeness and regularity theory which constitutes the basis 
of our analysis. Theorem 12.61 contains the statement of the novel regularity results 
established for the fluid component of the system. 

2.1 The PDE problem 

Let 0/ and £l s be the open smooth (three-dimensional) domains occupied by the fluid 
and the solid, respectively. The boundary of ft s represents the interface between the 
fluid and the solid, and is denoted by T s = dQ s ; Tf is the outer boundary of flf, 
that is Tf = dO, f \ d£l s . It is assumed that the motion of the solid is entirely due to 
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infinitesimal displacements, and hence that the interface F s is stationary. The PDE 
system under investigation comprises a linearized Navier-Stokes system for the fluid 
velocity u in O/, and an elastic equation for the displacement w of the solid in 
The coupling takes place on the interface T s ; the case of small, but rapid oscillations 
of the elastic body is considered, which yields the interface condition u = w t on T s in 
lieu of the no-slip boundary condition u = 0. (For a detailed discussion of modeling 
issues see, e.g. [H].) Thus, if p denotes the fluid pressure, the pair (u,w,p) satisfies 
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In the above PDE model e and <r denote the strain tensor and the elastic stress tensor, 
respectively, that are 

e y'M = 2 (^T + ^r) ' °»J M = A e fcfc (a;) + 2/i e^- (a;) , (2.2) 
^ 4 fc=i 

where A, u > are the Lame constants and 5y is the Kronecker symbol. The (uncou- 
pled) linear system for the fluid velocity field, often referred to as Oseen equations, 
is the linearization of the Navier-Stokes equations about a given (time-independent) 
smooth vector valued function v such that divv = 0; see, e.g., [T§]. Thus, the operator 
C is defined as follows: 

Cu(x, t) = (v(x) ■ V)u(x, t) + (u(x, t) ■ V)v(x) . (2.3) 

2.2 Functional setting, well-posedness 

The energy space for the PDE problem (|2.1j) is 

Y = HxH 1 (Q s ) xL 2 (n s ), 

where H is defined as follows: 

H := {u £ Lziflf) : divu = , u • u\r f = 0} , 

where the divergence is meant in a distributional sense. Ph will denote the Leray's 
projection operator from L2(f2/) onto H; then, L = Ph£- In addition, V denotes the 
space defined as follows: 

V := {ue H l {Vt f ) : div u = , u\ Vf = 0} ; 
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The following distinct inner products will be used throughout the paper: 

(u,v)f-.= / uvdQf, (u,v) s := / uv dQ s , (u,v) := / uvdT 
Jn f JQs ./I.. 

The space V is topologized with respect to the inner product given by 

(u,v) ljf := / e(u)e(v)dn f ; 
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the corresponding (induced) norm | • |i t is equivalent to the usual H^^lf) norm, in 
view of Korn inequality and the Poincare inequality. The norm || • || n r (D) m the Sobolev 
space H r (D) will be denoted by | • \ r< u throughout. 

Remark 2.1. Note that all the Sobolev spaces H r related to u and w are actually 
(H r ) 3 : the exponent is omitted just for simplicity of notation. 

Let us recall from [6] the definition of weak solutions to the (uncontrolled) PDE 
system (|2.ip . 



Definition 2.2 (Weak solution). Let (uo,wq,Zq) G H x ^(flg) x L2^l s ) and T > 0. 
We say that a triple (u,w,Wt) G C([0, T],H x i? 1 ($l s ) x L2(Q S )) is a weak solution to 
the PDE system ([HE]) if 

• (u{-, 0), w(-, 0),w t (-,0)) = (u 0l w ,z ), 

• u£L 2 (0,T;V), 

• a{w) ■ v G L 2 {Q,T-H- l / 2 (T s )), ±w\r s = u\r s G L 2 (0, T; H l l 2 (T s )), and 

• i/ie following variational system holds a.e. in t G (0,T); 

f •% + ( e (^)> e (^))/ + 0) / - (^M • ^ 0) = o 

1 |(^t, 1>), + (a(w),e{tP)) s - (a(w) ■ v, V) = , 
/or a// £esi functions (j) £ V and tp G .K^fis). 



(2.4) 



Existence of weak (global) solutions to the PDE problem (|2.ip follows by a more 
general — because pertaining to the actual nonlinear model, which comprises the Navier- 
Stokes system — result established in |6J. 

Remark 2.3. It is important to emphasise that the regularity of boundary traces 
of the elastic component required by Definition 12.21 does not follow from the interior 
regularity (u, w, wt) G C([0, T], H x H 1 (^ s ) x L2(Q S )) via standard trace theory. This 
is an independent regularity result, valid for any solution corresponding to initial data 
in the finite energy space; see the statement of Theorem 12.41 below. 
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Theorem 2.4 (Existence of weak solutions, [6]). For any initial datum (uq,wq, Zq) G Y 
and any T > 0, there exists a weak solution (u,w,wt) to system (|2,ip . along with the 
following trace regularity for the elastic component of the system: 

Vw\ Ts eL 2 (0,T;H-^ 2 (T s )), ^w\ ra = wt\ Ts € L 2 (0, T; i/ 1/2 (r s )) . 

It is important to recall from [6] that by properly defining the overall dynamics 
operator A, the PDE system (|2.ip admits as well a semigroup representation. Prelim- 
inarly, one needs to introduce the fluid dynamic operator A : V — » V, that is 

(Au,<f>) =-(e(«), £(</»)) V^ey, (2.5) 

and the (Neumann) map N : L2(r s ) — » H defined as follows: 

N-.g^h^(e(h),e(4>)) = {g,4>) V0 G F . (2.6) 

Then, by using the operators introduced above, the initial/boundary value problem 
(|2.1j) becomes the following abstract differential system in the variable y := (u,w,wt) 

y' = Ay, y(0) = i/o, (2.7) 

where yo := ( u 0; ^o> z o) £ ^ an d «4 : 2?(^4) C Y — » Y" is defined by 

I A-L ANo-(-) -v \ 
^l:= J J , (2.8) 

\ diver / 

with domain 

V(A) := {y = (u,w,z) eY : u G V, A(n + iVcr(w) ■ u) - Lu & H, z G i? 1 ^), 

div<7(«>) G L 2 (S7 S ) ; z|r s = u|r s } • (2.9) 

It was shown in [6, Proposition 3.1] that the operator A : T> (A) CY^Y defined by 
(|2.8jl - (|2.9p is the generator of a Co-semigroup e At of contractions on Y. Consequently, 
any semigroup solution y := (u, w, wt) belongs to C([0, T], Y) when y$ = (uq, wq, Zq) G 
Y, while it belongs to C([0, T],T>(A)) if yo G T>(A). A major achievement of the study 
performed in [6] is that it gives, in addition, that any semigroup solution turns out to 
be as well a weak solution of system (I2.ip . In particular, it satisfies the variational 
formulation (|2.4p . along with the intrinsic boundary regularity of the stresses on the 
interface. 

A different semigroup set-up, which does not require — at the level of (semigroup) 
generalized solutions — the present boundary regularity of the normal stresses, has 
been developed in [3]; see also [HE]. 

Remark 2.5. We finally recall from [2T| that the operator A defined by (|2.5h may be 
considered as acting on H with domain V{A) := {u G V : |(e(n), e(0))| < C|0|#}. 
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Readily A : T>(A) C H — > H is a self-adjoint, negative operator and therefore it 
generates an analytic semigroup e At on H. Then, the fractional powers of —A are well 
defined. A well known result relating (— A) a and the semigroup generated by A is the 
pointwise (in time) estimate, true for any < a < 1: 

\\{-ATe M \\ c[H) ~^, 0<t<T. (2.10) 

(To simplify the notation, we shall write A a , rather than (— A) a , throughout; similarly, 
later A a will mean (—A) a .) 

2.3 New results: regularity of boundary traces for the fluid compo- 
nent 

The main contribution of this work is a set of regularity results for the traces of the 
fluid velocity field u (and Ut) on the interface between the fluid and the solid regions, 
collected in Theorem 12.61 below. While these regularity estimates are of intrinsic inter- 
est, they additionally play a major role in the study of the quadratic optimal control 
problems associated with the PDE system under investigation, in the presence of con- 
trol actions on the interface. This fact will be briefly discussed in Section HI 

Theorem 2.6 (Main result: boundary regularity for the fluid component). Let y(t) = 
(u(t),w(t),wt(t)) be the solution to the coupled PDE system ()2.ip corresponding to an 
initial state yo = (uq,wq, zq). The following assertions pertaining to the regularity of 
boundary traces for the fluid velocity field u ( and ut ) are valid. 

1. If yo G Y, then 

u\ Ts EL^ h (0,T;L 2 (T s ))+L p (p,T;L 2 (T s )) Vp > 1 , (2.11) 

for arbitrarily small h > 0. More precisely, u admits a decomposition 

u(t) = Ul (t) + u 2 (t) , (2.12) 

where u\ satisfies a pointwise (in time) "singular estimate", namely there exists 
a positive constant Ct such that 

\\Mt)L 2 (r s ) < ^shoh Vy G Y , Vi € (0,T] , (2.13) 
with arbitrarily small 5 > 0, while 

u 2 \r s eL p {0,T;L 2 {T s )) Vp > 1 . (2.14) 

2. If yo G V{A e ) for some e G (0, \), then the regularity of u 2 in (|2.14p improves 
as follows: 

u 2 \r s GC([0,T},L 2 (T s )). (2.15) 
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3. Let now yo £ T>(A 1 9 ), with 9 £ (0, \). Then, there exists q £ (1,2) such that 

ut\r, EL q (0,T;L 2 (T s )) (2.16) 
continuously with respect to yo; namely, there exists a constant Ct such that 

IKIIz 9 (o,T;Z 2 (r s )) < C T WvoWviA 1 - 9 ) • ( 2 - 17 ) 
The exponent q will depend on 6: more precisely, given 6 £ (0, 4), we have 

2.4 Regularity of the boundary traces for the elastic component 

We conclude this section by introducing a key regularity result for the normal stresses 
of the elastic component of the system and a trace estimate pertaining to the fluid 
velocity field (Lemma 12.71 and Lemma 12.81 below, respectively), which are central to 
the proof of our main result. We fashion the proof of the assertions in Lemma 12.71 in 
some details, following the arguments formerly used in [2H Lemma 5.2]. Lemma 12.81 
will follow accordingly. 

The first result stated below is a slight variant of the one given in [211 Lemma 
5.2]. More precisely, while the regularity established in |21|, Lemma 5.2] pertains to 
the uncoupled wave equation and is expressed in terms of Sobolev norms, the result 
established in Lemma 12.71 concerns the regularity of the elastic component of the 
coupled system driven by the semigroup e" 4 *. Thus, the regularity of initial data is 
naturally measured by means of norms in intermediate spaces between the domain 
of the semigroup generator T>{A) and the state space Y; see (|2.2ip below. (Recall 
from Remark 12.51 that for simplicity of notation we write A a instead of (—A) a .) It 
is this latter formulation, more appropriate in the present context, to require some 
adjustment of the arguments used for the proof of Lemma 5.2 in [21]. The technical 
details are the same as the ones in [21], if one takes into account that the solutions 
considered are generated by the semigroup flow. However, for the reader's convenience 
the proof of Lemma 12.71 follows below. 



Lemma 2.7. Let (u,w,wt) be the solution to the PDE problem (12. ID corresponding 
to initial data yo = (uq,wq,zo) £ V(A a ), < a < 1, and let f := u\r 3 so that 
f £ £2(0, T; i? 1//2 (r s )). Then, the w component, which solves the initial/boundary 
value problem 

wu — divcr(tt;) =0 in Q s 

£ t w\ Ts =f on (2.19) 

w(-,Q) = w , w t (-,0) = zo in tt s , 

can be decomposed as w = w\ + W2, where a(w\) ■ v £ C([0, T], i/~ -1 / 2 (r s )) ; while 
a{w2) • v £ ^2(0, T; L2(r s )). //, in addition, f £ H a (L s ), then cr(w2) ■ v £ H a (T, s ). 
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Moreover, the following estimates hold true. 



IkOl) ■ I/ llc7([0,T],iJ-- 1 /a(r s )) - C l{\ w Q\l,Sls + 1*0 1 0,0, + ll/lli 2 (0,r;HV2(r s )) ( 2 - 20 ) 

\WM ■ Hl!r«(E.) ^ C 2 (\\A a y \\ H + \\f\\ H «p.)) (2-21) 

Proof. The first part of the Lemma follows directly from [6, Theorem 3.3]. Similarly, 
the regularity inequality (|2.20p is implied by [21\ Lemma 5.2]. Indeed, this estimate 
involves the energy level regularity of the data (wq,Zo), and therefore it is a simple 
consequence of a more general result — valid for any open loop problem like f)2. 19|) — 
established in |2H Lemma 5.2]. The inequality in (|2.2ip involves instead initial data 
whose regularity is 'above' the energy level, and it is thus expressed in terms of the 
domains of fractional powers of the generator A (rather than fractional Sobolev spaces 
as in [21]). 

In order to establish (12.21H we appeal, as in [21], to interpolation methods. Ac- 
cordingly, it will suffice to prove (|2.21j) for a = and a = 1. The regularity result 
corresponding to a = has been shown in [6, Theorem 3.3]; thus, it remains to an- 
alyze the case a = 1. Note that (|2.21l) with a = 1 provides IT (£ s )-regularity of the 
component a{ui2) in terms of ||.4yo||.ff + ||/||j3"i(s s )- As in [21] and [6], we exploit the 
fact that W2 = x w i where x is supported in the microlocal hyperbolic sector; see [U 
Section 5] . This means that the tangential boundary regularity is implied by the time 
regularity alone. In other words, if 1 (E s )-regularity is equivalent to H 1 ^, T; L2(T S )) 
for functions supported in the hyperbolic sector; see [61 Section 5] for more details. 

In view of the above remarks, we need to show the following estimate: 

IffW T ^ < C(\\Ay \\ H + l|/IUx (Es) ) . (2.22) 

In turn, this inequality expresses the so called "hidden regularity" for the wave equa- 
tion, written in the variable w = wt, where u)(0) = ZQ,wt{0) = div<r(u;o), w\r 3 = /• 
Namely, w satisfies 

wtt — div a(w) = 

y w(-,0) = z , w t (-,Q) = divo-(wo) . 

Owing to Theorem 2.1 in [22], we obtain 

Mw) ■ u\\ L2 p 3 ) < C(||/|Ux (Es) + \\w(;0)\\hHU) + 0)11^(0)) , (2.23) 

subject to the compatibility condition w(-,0)\r 3 = /(-,0). The conditions needed for 
applicability of the said regularity are 

U G H 1 ^) , dwa(w ) G L 2 (n s ) , / G H 1 (E 8 ) , 

< (2.24) 

[with the compatibility condition £o(vlr s = /(-,0). 
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Then, ([233]) reduces to: 

■ f|U a (E.) < c (\\f\\m(E 3 ) + l|divo-(™ )||L 2 (n s ) + ||*b||ffi(n,)) • (2.25) 

Notice now that the conditions (|2.24p are satisfied by any yo = (uo,wq,zo) G D(A); 
indeed, from the very definition (|2.9p (of the domain of the generator A) it follows 

z GH l (Q), diva(w )£L 2 (n s ), z \ Tb = u(-, 0)|r s =/(-,0), 

where the last equality follows from the identification / = u|p s . This yields 

Hw) ■ u\\ L2{Ss) < C (||/|bi (Ss) + \\Av4h) , (2.26) 

which in turn implies (|2.22p . after restricting w to x w = w 2 (we use here the same 
arguments as in [6]: u>2 = X w satisfies the wave equation with microlocal support 
in the hyperbolic sector). Interpolation with the inequality corresponding to a = 
finally establishes ([23T]) . □ 

The sharp results of Lemma 12.71 yield the following improved boundary regularity 
for the fluid velocity field. 

Lemma 2.8. Consider the PDE problem (|2.ip , and take initial data = (uq, wq, zq) 6 
T>(A a ), < a < -7. Then, for any T < oo we have u\r s € H a CE s ) and the following 
estimate holds true: 

IMIff»(E.) <C\\A a y \\Y. (2.27) 

The proof of this Lemma is identical with the one given in |21|, Lemma 5.3], after 
taking into consideration the estimates derived in Lemma 12.71 



3 Proof of the main result 

This section is devoted to the proof of Theorem 12.61 As we shall see, the regularity of 
boundary traces for the fluid component of the system will follow combining the excep- 
tional boundary regularity of the hyperbolic component of the system (Lemma 12 .7p . 
the smoothing properties of the (fluid) analytic semigroup, and the theory of interpo- 
lation spaces. 

To begin with, let us recall from [21] the basic regularity of the Neumann map, 
which will used repeatedly in the proof of our main result. 

Lemma 3.1 ([21]). For the map N : L2(T S ) — > H defined in (|2.6p the following 
regularity results holds. 

(i) N G C(L 2 (T s ),V(A 3 / 4 - s )) nC{H- x l 2 {T s ),V{A 1 l 2 )) for any 5, < 5 < |. 

(ii) One has N*Au = — u\r s , u G V, where the adjoint of N is computed with respect 
to the L2-topology. 
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This Lemma is well known in the context of classical elliptic equations, where iV 
is the associated Green (Neumann) map. The same result has been proved for the 
abstract operators (defined by (|2.5p and (|2.6p ) which arise in the PDE problem under 
investigation; see |214 Proposition 4.3]. For the reader's convenience we recall few 
details. 

The composition N*A (N* being the L2-adjoint of N) is computed from 

(N*Au,g) = (Au,Ng) f = -(u\ Ts ,g) for all u G V(A), g G L 2 (T S ). 

Since u|r s is bounded on V = T>(A 1 ^ 2 ), the above identity can be extended by density 
to all of V, even more, to V{A l / i+s ). This gives the identification in part (ii) of the 
Lemma. 

Remark 3.2. In view of part (i) of Lemma 13.14 we can define the operator AN as 
acting from L 2 (T S ) into [P(^4 1 / 4+5 )] / . On the other hand, the analyticity of e At yields 
the standard bound ||^e^*||£(m < Ct~ s , < t < 1. This estimate, combined with 
the invariance of the semigroup e At on the domains of fractional powers of A, allows to 
define the convolution e A ^~ s ^ ANg(s)ds with values (for each t) in T>(A 3 / A ~ S ) C V, 
for any g G L oo (0,T;L 2 (T s )). Thus, —N*A /* e^ 1 '^ ANg(s)ds can be identified— on 
the strength of Lemma I3.H — with the trace on the boundary of e A ^~ s ^ AN g(s)ds. 
In the remainder of this section we will use freely the above identifications. 

Proof of Theorem \2.6l Our starting point is the equation satisfied by u(-), namely 
ut = (A — L)u + ANa(w) ■ v, whose evolution/integral equivalent is 

u(t) = e M u - f e A{t ~ s) Lu(s) ds + I e A{t ~ s) ANa{w{s)) -uds. (3.1) 
Jo Jo 

We claim that the expected decomposition (|2.12p . along with the various regularity 
properties listed in the statement of Theorem 12.61 is the one with 



ui(t) := e At u , u 2 (t) := - f e A{t ~ s) 'Lu(s) ds + f e A(t ~ s) ANa{w{s)) -uds. (3.2) 

Jo Jo 

1. Since by Lemma 13. II we know that n|r s = —N*Au, the splitting (13.21) yields as well 



Ul | rj = -N*Aui(t) = -N*Ae At u , 

u 2 \ Ts = -N*Au 2 {t) = N*A [ e A{t ' s) Lu{s) ds - N* A [ e A(t - s) ANa(w(s)) -uds. 

Jo Jo 



Ui(t) U 2 (t) 

(3.3) 

respectively. The term u\ yields a pretty straightforward regularity estimate. Indeed, 
combining the (analytic) estimates (|2.10p with the regularity of the Neumann map 
recalled in Lemma |3.1| (ii), we obtain 

\\ui(t)\\L 2 (T s ) = \\N*Ae At u \\ L2{Ts) = \\N*A 3 /*- s AV*+ s e At u \\ L2{Ts) 

(3.4) 

< \\N*A 3 / 4 - 5 \\ \\A^+ s e M u \\ H < Ct-^WuoW < dt-^WyoWy , 
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1. e. the "singular" estimate f)2 . 13|) holds true. In particular, 

uijr s G L4_/j(0, T; L 2 (T S )) with arbitrarily small /i > 0. (3-5) 

Thus, to establish ()2.1ip . it remains to demonstrate that the boundary trace of the 
second component u 2 on T s is well defined and satisfies (|2.14j> . Taking account of 
(|3.3fl . since the regularity of U 2 has been dealt with in [12], we will be done once we 
show 

U 1 eL p {0,T;L 2 {T s )) Vp > 1 . (3.6) 

We consider the convolution integral z(t) = J * e A ^~ s ^ Lu{s) ds, and recall from [61 The- 
orem 3.2] the additional property u G 1/2(0, T; V), which implies Cu G £2(0, T; £2^/)) 
and Lu = PhC-u G £2(0, T;H). Thus, by parabolic regularity theory we deduce 
z G L 2 (0,T;V(A)) n i? 1 (0, T; if); therefore, by interpolation, 

2 € H a (0,T;T>(A 1 ~ a )) Va G (0, 1) . 

In particular, we may take a = | — 5 to find 

J8 G ^3/4-5(0^.^^1/4+^^ 

and since 

?7i(t) = [7V*yl 3 / 4_5 ] A 1 / 4+(J z(i) , 

we obtain 

Ux G ff 3/4-5 (0j T . L2(Fs)) c c([0) T ], L 2 (r s )) , (3.7) 

where the inclusion follows by Sobolev embedding theorems in one dimension. Thus, 
the validity of (|3.6p is ensured by (13. 7p . 

The analysis of the summand U 2 (t) critically relies on the sharp regularity theory 
for the elastic component of the system, as given by the first part of Lemma 12.71 
Specifically, on the basis of 

a(w) ■ v G C([0, T], H'^iTs)) + L 2 (0, T; L 2 (T S )) , 

it was shown in [12J that 

U 2 eC([0,T],L 2 (T s )) +L p (0,T;L 2 (r s )) C L p (0,T;L 2 (T s )) Vp > 1 . (3.8) 

Consequently, combining (|3.8p with (|3.6p we obtain (|2.14j) . as required. This concludes 
the proof of all assertions in item 1. 

2. Let now yo G T>(A € ), e > 0. Aiming to show the validity of (|2.15p . we recall 
u 2 \r s = —N*Au 2 = U\ — U 2 , with U{ defined in ()3.3p . We have seen in step 1 that 
Ui G C([0,T], L 2 (T S )) when 7/0 just belongs to Y; the same holds true, a fortiori, in 
the present case. 

As for U 2 , its regularity analysis when yo G V(A e ) has been pursued in [T2]. We 
just note that since 

yo e ^(^ e ) , 
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according to Lemma [2T81 u G H £ (T, S ), provided that e < \. Therefore, we may apply 
the second part of Lemma 12.71 (/ = w|r s ), thereby obtaining 



(3.9) 



Owing to (|3.9p . a decomposition of U 2 is necessitated, say U 2 = U21 + £^22- Careful 
calculations made in |T2] showed that now both terms U 2 i belong to C([0, T], L 2 (T S )), 
i = 1,2, which combined with (|3.7p finally yields the conclusion (|2,15p . 

3. We turn our attention to the time derivative Uf. In order to prove (|2.16p . we 
compute from (|3.ip (in the dual space [X>(A)]') 



u t (t) := Vl (t) + v 2 (t) = Ae At u + v 2 (t) , 
where, initially, yo G Y and v 2 reads as 

v 2 (t) = - [Lu(t) + A I e^-^Luis) ds 



(3.10) 



+ 



A 



jA(*-«) ANa(w(s)) ■ v ds + ANa(w(t)) ■ v =: -V x {t) + V 2 (t) . 



As V\ is nothing but Zt (z is defined in step 1), and according to parabolic regularity 
z G H l (0,T; H), we immediately obtain 

V 1 eL 2 (0,T;H)cL 2 (0,T;L 2 (n f )). (3.11) 

The more challenging analysis of V 2 yields, utilizing once again Lemma 12.71 

V 2 G L 2 (0,T; [Z>(AV2+*)]') = L 2 (0,T; [H 1+2S ^ (fy)]') , (3.12) 

for any < Si < ^; see |12j for all details. Combining (|3.1ip - (|3.12p leads us to the 
estimate 

y GY ^v 2 a 2 (0,T;[ff 1+Ml (%)]'), < <5i < J . (3.13) 
When yo 6 "D(A), we rewrite U2 in a different fashion, namely 

U2(t) = - \ e At Lu + I e A ^Lu s (s) ds 



+ 



e Mt-°) ANa(w s (s)) ■ v ds + Ae At Na(w(0)) ■ u] =: - Zi(t) + Z 2 (t) . (3.14) 



According to well-posedness theory in [6j Theorem 3.2] yo £ ^G^) implies now y^ = 
(ut,wt,wtt) G C([0, T];Y), with the additional regularity u t E L2(0, T;V). Conse- 
quently, Lut G L 2 (0,T; H) and 



e A (- s 'L Us (s) G L 2 (0, T; V{A)) C L 2 (0, T; F 2 (fy)) . (3.15) 
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On the other hand, with yo G *D(A) we have uq G V so that Luq G and in particular 

e A Lu G L 2 (0, T; ^(A 1 / 2 - 5 )) = L 2 (0, T; i^ 2 * 5 ^/)) , 
which is weaker than f)3. 15|) : consequently, 

Z 1 €L 2 (0,T;H 1 - 25 (n f )), < 5 < i . (3.16) 

Consider next the term Z 2 , and notice that yo € ^(-4-) implies, in particular, 
(ujo, Zo) ^ H 2 (Q S ) x By standard semigroup arguments and according to 

Lemma 12.71 we know that 

a(w t ) ■ v G C([0, T\,H-V 2 (r a )) + L 2 (0, T; L 2 (r s )) . 

In light of the above result, it was shown in [12J that the convolution integral 

Z 21 := f e A{ - s) ANa(w s (s))-vds 
Jo 

satisfies 

Z 2 i G C([0, T], D^ 1 / 2 - 5 )) + L 4 (0, T; L>(A 1//2 ~ 5 )) , < 5 < ~ , 
which gives 

Z 2 i GL 4 (0,T;Z)(A 1 / 2 - 5 )), 0<5<i. (3.17) 

On the other hand, still with y G D(A), one has just a(w)-v G C([0,T],.ff _1 /2(r g )). 
Rewrite the summand Z 22 (t) := Ae^iVc^u^O)) • v in (|3.14p as 

Z 22 (t) = A l / 2 e M [A l / 2 Na{w{<S)) ■ u] , 

and use once more Lemma [3. 1\ along with the singular estimates pertaining to analytic 
semigroups, to find 

Z22 G L q (0, T; D(A l/2 - 6 )) , provided that q(l - 5) < 1. (3.18) 

for any 5 < Recalling v 2 = — Z\ + Z 2 \ + Z 22 , with the requirements (|3.16j) . <|3. 1T|> 
and (|3.18p . we finally obtain 

y 6 P(.4) f 2 G ^^^^(i^-fe)) = L q2 {^T-H l ~ 2 ^(n f )) (3.19) 

for any < 5 2 < \, where q 2 G (1,2) depends on <5 2 ; more precisely, 

92 < t— V • ( 3 - 2 °) 
1 — 2 

We proceed now as in the conclusion of step (Hi) in |12} Proof of Theorem 2.9]. 
Combining fj3. 19j) with (|3.13p . gives by interpolation 

y G D(A 1 ~ e ) v 2 G L q2 (0, T; W) , (3.21) 
14 



where g 2 is as in ()3,20p and W is the interpolation space 

if 

s = (1 - 0)(1 - 25 2 ) - 0(1 + 2<*i) = 1 — 2tf 3 - 20(1 + <5i - 5 2 ) > ; 

see |241 Chapter 1, Theorem 12.5]. Notice that by taking, for instance, 6% = 62 ='■ 5, 
one has s > 1/2 provided that 

9 + 8 < - . (3.22) 
4 

In this case t> 2 takes values in H s (£lf) with s > 1/2 and hence its trace on T s is well 
defined. Notice that, in view of the constraint (|3.22p . we need to require < 9 < |. 
Consequently, given any 9 such that < 9 < j, choosing, e.g., 5 = 1/4 — 9 in view of 
rf3Tm from (I33TD it follows 

y G Z^ 1 - 9 )) =^ N*Av 2 e L q2 (0,T;L 2 (T s )) Vg 2 < . (3.23) 

It remains to establish the regularity of the first summand N*Av\(t) = N* Ae At Au$ 
when y G D(A 1 ^ e ). In this case u G (iT 1 ^/), L 2 (£lf))g = iT 1_e, (0/), and it is not 
difficult to show 

y e DiA 1 - 9 ) => N*A Vl G L qi (0, T; L 2 (r s )) V 9l < - - J - - (3.24) 

where e is independent of 9 and can be taken arbitrarily small. (Notice that the 
Sobolev exponent q\ belongs to (1,2), as well.) Thus, in view of (13.231) and ([3. 24ft we 
finally obtain 

y e DiA 1 - 9 ) =^ u t \r s GL 9 (0,T;L 2 (r s )) y q < —L- , (3.25) 
which concludes the proof. □ 



4 Application to the quadratic optimal control problem 

In this section we briefly discuss the implications of the obtained trace regularity results 
in the study of the optimal control problem (with quadratic functionals) associated 
with the fluid-solid interaction model (|2.ip . Thus, we consider the initial/boundary 
value problem 



ut — div e(u) + Cu + Vp = 


in 


Qf 


= n f 


x (0,T) 


divu = 


in 


Qf 






wtt — div a{w) = 


in 


Qs 


= n s 


x (0,T) 


u = 


on 


£/ 




x (0,T) 


w t = u 


on 


z s 


■■=r s 


x (0,T) 


a{w) ■ v = e(u) ■ v — pv — g 


on 


s s 






u(0, •) = Uq 


in 


n f 






w(Q, •) = w , w t (0, •) = z 


in 


£i s , 
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which includes a function g acting as a control on the interface T s . The PDE problem 
(|4.ip readily corresponds to an abstract (linear) control system, that is 



\y' = Ay + Bg in [V(A*)Y (42) 
\y(0)=y eY 

where the (free dynamics) operator A is the one in (|2.8p . while the (control) operator 
B : L 2 (r s ) =:U^ [D{A)}' is given by 



B = . (4.3) 




Semigroup well-posedness of the uncontrolled problem corresponding to (|4,2|) still fol- 
lows from [6l Proposition 3.1]. The membership B 6 C(U, [D(A)]') is a consequence 
of AN being bounded from L 2 (T S ) into V (yet unbounded from L 2 (r s ) into -ff); see 
|21j for more details. 

We associate to (|4.2p a quadratic functional 

J(g)= I {\Ky(t)\l + \g(t)\l)dt (4.4) 
J o 

where (7 varies in £2(0, T; £7), and initially 7?- € C(Y,Z) (7Z is called observation op- 
erator, Z is the observations space). Then, the (optimal) quadratic control problem 
over a finite horizon is stated as follows. 

Problem (The optimal control problem). Given yo € Y , seek a control function 
g G L,2(0,T;U) which minimizes the cost functional (14. 4p . where y(t) = y(t;yo,g) is 
the solution to (|4.2p corresponding to g. 

As in the classical Linear Quadratic (LQ) problem in a finite dimensional context, 
solvability of the optimal control problem defined above does not just mean existence of 
a (unique) minimizer g* . We aim to obtain the feedback representation of the optimal 
control, defined for any 7/0 hi the state space Y ', along with well-posedness of the 
corresponding operator Riccati equations. On the other hand, it is well known that 
in the case of boundary control systems, such as the present (|4.2p . a major difficulty 
comes from the fact that the gain operator B*P(t) which occurrs in the feedback 
formula 

g*(t) = -B*P(t)y*(t) a.e. in[0,T], 

may not be defined, or may be unbounded. 

Let us recall that owing to the well-known representation of the Riccati operator 
P(t) (involving the evolution map <£(•, •) associated with the optimal dynamics) 

r T 

P(t)z = J e A * {s - t) R*R<Z>(s,t)zds, z£Y, 
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it is immediately seen that boundedness of the operator B*P(t) will crucially depend 
upon the regularity properties of the operator 

B*e AH : V(A*) ^U = L 2 (T S ) . 

Thus, a PDE intepretation of the composition B*e A * is called for. This is obtained — 
for a given control system (|4.2p — by computing B* and by noting that e A * l yQ solves 
the uncontrolled Cauchy problem y' = A*y, y(0) = yo. These calculations eventually 
yield the distinct trace operator corresponding to B*e A * t . (This connection appeared 
clear since the former studies of the LQ-problem for simple hyperbolic equations; see, 
e.g., [18J and [8].) In the present case we obtain, in particular, 

B*e A *-y = [N*A, 0, 0]y = N*Au = -u\ Ts , (4.5) 

where y = (u,w,wt) solves the initial/boundary value problem corresponding to the 
afore-said adjoint system, which is almost identical with the PDE problem (|2.1j) . This 
explains why a study of the boundary regularity of the fluid velocity field on the 
interface is necessitated. 

A situation where one can conclude that B*P(t) is bounded on the state space Y 
is when the pair (A, B) yields a 'singular' estimate for the operator e At B near t = 0; 
see |17j . [19] . and [20]. It was shown in [21] that the PDE problem (|4.ip does give rise 
to a singular estimate; however, this holds true provided the observation operator 1Z 
possesses a suitable smoothing property. 

In contrast, in light of the regularity results established in Theorem 1 2. 6\ we are able 
to show that the abstract control system (|4.2p corresponding to the initial/boundary 
value problem (|4.ip falls in the more general class of systems introduced in [2j. In- 
deed, all the regularity results proved in Theorem 12.61 — which yield relative assertions 
for a similar (dual) PDE problem — are the PDE counterpart of the control-theoretic 
properties (Hypotheses 2.1-2.2 therein) required to apply the theory of [2]. We do 
not discuss the proof, which follows the lines of [X2^ Proof of Theorem 2.6], but let 
us remark that while the regularity of the boundary traces of u on T s corresponds to 
the regularity of the operator B*e A * t , the obtained regularity for Ut in (|2.16p yields a 
sought-after and most challenging property for B*e A * t A* (that is condition (3.23) in 
[2], or (4.4) in [32]). We just note that the latter PDE interpretation follows from the 
following key observation: 

B*e A * t A*% = B*j t e AH {A* e ~ l yv) = N* Au t = -u t \ Ts , 

where (u, w, wt) solves the same (adjoint) PDE problem as before, yet with initial state 
£o = A* e ~ 1 yo € V{A* 1 ^ 9 ) (while time derivatives are understood distributionally) ; the 
reader is referred to [12 , Section 4.1] for more details. 

Thus, in order for all the hypotheses needed to apply [2J Theorem 2.3] to be fulfilled, 
we just need to assume the one which pertains to the observation operator 1Z; namely, 

11*11 G C(V(A e ), V(A* e )) (4.6) 
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for some e € (0, 1/4). Then, under (|4.6p all the conclusions of [2j Theorem 2.3] follow. 
In particular, one obtains the feedback synthesis of the optimal control for all y$ in the 
finite energy space Y, with a gain operator which is densely defined, that is enough to 
get well-posedness of the corresponding differential Riccati equations. 

Remark 4.1. We finally note that while a case of particular interest is 1Z = I, 
the requirement (|4.6p does not 'force' smoothing effects of the observation operator. 
Indeed, in the present case T>(A* e ) coincides with V{A e ) for e sufficiently small, and 
hence (I4.6P is satisfied by any operator which 'maintains' regularity, such as the identity 
1Z = I; see \12\ Remark 2.7] for a more detailed explanation. This answers the final 
question raised in [211 Remark 6.1]. 
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